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Numerous studies have demonstrated the potential for simple fluid plus particle systems to pro- 
duce complicated dynamical behavior. In this work, we study a horizontal rotating drum filled with 
pure glycerol and three large, heavy spheres. The rotation of the drum causes the spheres to cascade 
and tumble and thus interact with each other. We find several different behaviors of the spheres 
depending on the drum rotation rate. Simpler states include the spheres remaining well separated, 
or states where two or all three of the spheres come together and cascade together. We also see 
two more complex states, where two or three of the spheres move chaotically. We characterize these 
chaotic states and find that in many respects they are quite unpredictable. This experiment serves 
as a simple model system to demonstrate chaotic behavior in fluid dynamical systems. 

PACS numbers: 47.32.Ef, 47.52. +j, 47.80.Jk 



I. INTRODUCTION AND PRIOR WORK 

Prior experiments and simulations have shown that 
particles interacting with fluids can behave chaotically. 
Simple fluid flows with nonspherical particles 1, 2], 
stirred particle-filled fluids [||, and various other hydro- 
dynamic interactions Q can all produce chaotic behavior. 
In this manuscript, we look at chaotic behaviors involv- 
ing particles in a horizontal drum which cascade under 
the influence of a constant external force (i.e. gravity), 
as first studied by Mullin et al. [5j]. 

In the study of classical Newtonian chaotic systems, 
there are several heavily studied systems which serve as 
examples to which other systems can be compared. Sim- 
ple experiments such as a forced pendulum (fj are easily 
constructed laboratory experiments which can be stud- 
ied in depth to better understand low dimensional chaotic 
behavior. Fluid systems, while containing many more de- 
grees of freedom, also exhibit several geometrically simple 
experiments which can produce interesting results, giv- 
ing us insight jnto the character of chaotic fluid dynam- 
ical systems (7|. Relevant to this manuscript, consider 
non-Brownian particles sedimenting under the influence 
of gravity. When many particles sediment, interesting 
swirling motions are seen 0. When only a few parti- 
cles sediment interesting effects are still seen, including 
chaotic behavior [tj [h| • Specifically, if three particles are 
released close to each other, two of the particles eventu- 
ally pair up and travel together, while the third particle 
moves on its own isolated trajectory. Intriguingly, this 
pairing is sensitive to initial conditions. A slight change 
in the initial particle positions can lead to a switch of 
which particle is isolated. 

These simple sedimentation experiments have only 
transient chaotic behavior which is difficult to study ex- 
perimentally. By considering motion in a rotating drum, 
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FIG. 1: Sketch of rotating drum filled with viscous fluid and 
three beads. These beads are dragged upwards by the front 
wall of the drum until gravity pulls them away from the wall 
and they fall through the fluid. (Based on Ref. 



particles can interact repeatedly, leading to longer- lasting 
and potentially persistent chaotic behaviors. In a rotat- 
ing drum partially filled with a liq uid and particles, seg- 
regation eventually occurs [ill Il2| . In a purely granular 
system (no added liquid), similar segregation is also seen 
0, [3, [H, [H, [l?], EH . A simpler system was devised by 
Mullin et al., in analogy with the three-particle sedimen- 
tation simulations of Ref. [n|. They studied a hollow 
cylinder, filled with glycerine, oriented horizontally, and 
rotated at various angular speeds u> @ . Three large non- 
Brownian spheres are placed in this drum. This system 
serves as a version of a simple sedimentation experiment, 
with the constraint that the rotation of the drum forces 
the three beads to remain close together, continually in- 
teracting. Figure [T] shows a sketch of this experiment. 
When the drum rotates, these beads can undergo a pe- 
riodic cascade in the vertical direction. As they cascade, 
the beads experience long-range interactions due to the 
fluid. Mullin et al. found several different behaviors de- 
pending on uj. Simple behaviors included states where 
the three particles remain well-separated and cascade in- 
dependently of each other. More complex behaviors were 
observed, such as a state described as "chaotic" where the 
particles move back and forth horizontally along the tube 
as they continue to tumble vertically. These horizontal 
motions were slow, erratic, and unpredictable. In some 
cases two particles even momentarily collide before with- 
drawing. No chaotic behavior was observed when there 
were only one or two beads. 
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Motivated by Ref. [5J, we study a similar experimen- 
tal system of a rotating drum with three spheres. The 
benefit of this experiment is that it can be stably run for 
very long periods of time, thus allowing us to study the 
motions of the spheres over a long period of time and 
describe the chaotic states in great detail (which was not 
done in Ref. []|). We find the states seen in Ref. [j| as 
well as three new states. Additionally, we characterize 
the chaotic states, finding that different chaotic states at 
different rotation rates u> often have widely different be- 
haviors, despite a superficial similarity. We thus demon- 
strate and characterize a simple system which possesses 
a wide range of nontrivial behaviors. 



II. EXPERIMENTAL METHODS 
A. The Drum 

The experimental apparatus is a / c = 25 cm long 
horizontally oriented sealed drum, with an inner radius 
r c = 5.72 cm. These dimensions are similar to those 
used in Ref. [|[ (l c = 25 cm and r c = 5.9 cm). The main 
body of the drum is constructed from a section of acrylic 
glass pipe with a wall thickness of 1 cm. To each end of 
this drum, waterproof threaded aluminum caps are fitted. 
Each cap has a shaft attached via an adjustable mount- 
ing, so that the shaft can be carefully centered within 
the cap. These shafts attach to a stand via two bearings, 
which allow the drum to rotate freely about the hori- 
zontal axis. The base of the drum stand contains four 
adjustment screws for leveling the apparatus. 

On one of the drum shafts, a pulley is mounted, and 
connected via a belt to a Dayton 1/2 HP 3-phase A/C 
motor driven by a Fuji AF-300 controller. This allows the 
drum to rotate on its axis at a variable rotation rate ui 
from 5 — 13 rad/s. Due to the belt drive connection, the 
actual motor rotation rate is potentially different from 
that displayed on the control box, so the rotation rate 
of the drum is measured independently using a Pasco 
PS-2120 rotary motion sensor connected to a computer. 
Measurements taken over four hours show the rotation 
rate to be stable to within 1%. As all experiments are 
started from rest, we also measure the time for the drum 
to spin up to within 1% of its final velocity, and find 
typical spin- up times to be on the order of tens of seconds. 

The drum is filled with 99.5% pure glycerol from 
Sigma- Aldrich, and three 440c stainless steel ball bear- 
ings purchased from Winstead Precision Ball Company, 
each of which has a diameter of 2r^ = 1.59 cm, density 
Pb = 7.65 x 10 3 kg/m 3 , and mass of 16.1 grams. When 
immersed in glycerol (pt = 1.26 x 10 3 kg/m 2 ), the beads 
have an apparent buoyant weight of 0.131 N each. 

To maintain constant viscosity, we control the temper- 
ature of the fluid within the drum (l9j |. To do so, we 
immerse the drum in a tank of water, in which a cop- 
per heat exchanger has been placed. We connect this 
heat exchanger to a Thermo NESLAB RTE-7 digital re- 



frigerated/heated bath, which is maintained at 25° C 
in all experiments. The rotation of the drum provides 
sufficient mixing to allow the water within the tank to 
be maintained at 24.7±0.1° C, as measured by a dig- 
ital thermometer, independently for each experiment. 
This results in a measured kinematic viscosity v = 7.69 
Stokes (as compared with 9.36 Stokes for the fluid used 
in Ref. [|). 

B. Data Collection 

To image and track the particles, we use a Pixelink 
PL-B741F 1.3 megapixel firewire monochromatic cam- 
era connected to a personal computer running Windows 
XP. A Monarch Instruments Nova-Strobe DAX is used 
to light the particles from behind, which minimizes the 
motion blur. An opaque screen surrounds the entire ap- 
paratus to block ambient light. Compressed movies, last- 
ing up to six hours in duration, were captured using a 
custom application written in CH — V and analyzed with a 
particle tracking algorithm implemented in Matlab. We 
can locate particle positions with a resolution of ±0.25r;,, 
limited by slight optical distortions. 

For each experiment, we initialize the particle positions 
by setting the drum rotation rate lu to that of a known 
chaotic state, and stopping the drum when the three 
beads are distributed equidistant from one another, with 
approximately 3.5 particle diameters spacing between the 
beads. The exact rotation rate is not important, as it is 
simply used as a tool to position the particles. Once the 
drum is stopped and enough time allowed to elapse for 
any fluid motion to cease (~ 10 min), we start the ro- 
tation of the motor at the desired rotation rate w, and 
immediately begin recording video. Videos are streamed 
directly to the PC hard drive into a compressed Xvid 
MPEG-4 AVI file. Video files are then post-processed 
using Matlab. Note that all experiments discussed follow 
this protocol, starting the drum from rest, setting the 
speed on the motor for the desired u>, and then turning 
the motor on. In particular, we do not examine hysteretic 
effects (although we speculate that there likely are some 
hysteretic effects, as discussed below). 

C. Sphere Motion and Nondimensional Numbers 

When the drum first begins to rotate, there is a tran- 
sient state where the fluid is not yet equilibrated to the 
new rotation rate. This can be estimated by the Ekman 
pumping time, 

TE = -^=, (1) 

based on the cylinder length l c , viscosity v, and final ro- 
tation rate uj [20(. For our experiment, with u> ~ 5 — 13 
rad/s, we have te ~ 2.4 — 3.7 s. This is faster than the 
time needed for the motor to reach full speed (tens of 
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seconds as noted above), and so the motor is the limiting 
factor in reaching the steady state. This implies that the 
behavior seen in our measurements, lasting several hours 
in duration, is not dependent on effects from the initial 
spin-up of the drum fluid. Observations of tracer mate- 
rial within the fluid confirm the short spin-up time. In 
Sec IIII D[ we will show that typical time scales for hori- 
zontal motion of the spheres (x direction) are 0(100 s), 
two orders of magnitude removed from the Ekman time 
rj ~ 3 s. The particle turnover time (the time for a 
single cascade to occur) is of order ~ 1 s. 
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FIG. 2: (Color online) Cascade speeds v up and Vdown plotted 
versus to. The solid line represents v — u>r c , showing that wr c 
is an approximate upper bound for the cascade speed. 



If we consider the ratio of inertial forces for the parti- 
cles compared to the fluid, we can typify the relative con- 
tributions by comparing the density ratio Pb/pf ~ 6.2, 
which implies that fluid inertia will have less influence 
upon the bead trajectories than the Reynolds number 
might otherwise imply. Due to their large inertia rela- 
tive to the fluid, the beads will not tend to follow fluid 
streamlines exactly, but rather interact via drag. 

Finally, we consider the Galilei number to quantify the 
ratio of gravitational to viscous forces. The force due to 
gravity is F g = (4/3)^1 Apg, and the viscous drag is 
defined as the Stokes drag F v = QnvrbU . Using these we 
define the Galilei number as [2l| 

with Ga = 3.9 for oj = 5 rad/s and Ga = 1.3 for w = 15 
rad/s, using the typical velocity scale as U — uir c . This 
shows that the influence due to gravity is always compa- 
rable to that due to viscous drag, which is not surprising. 
At lower rotation rates, the force due to gravity is pro- 
portionally larger, meaning the particles will not be lifted 
up as high in the cylinder (y direction) before falling back 
down; this is indeed what we observe. Likewise at higher 
rotation rates, the forces due to gravity and viscous drag 
are equal in magnitude, and the cascading motion carries 
the spheres further upward in y. 



When we are within the cascade regime, the beads re- 
peatedly cascade in the y — z plane, being dragged up the 
front wall of the cylinder, then falling away from it. We 
measure the average speed of the spheres as they rise and 
fall, and plot these speeds as a function of uu in Fig. [2] 
The solid line drawn at v = u;r c is an approximate upper 
limit to the characteristic speed U of the particles. Us- 
ing this, we can compute the particle Reynolds number 
Re = UL/v, which determines the contribution of fluid 
inertial effects relative to viscous drag. Using the bead 
radius to set the characteristic length scale L = r^, the 
Reynolds number is 



Re = 



ujr c r b 



i' 



(2) 



Noting that, for this experiment, r c , r?,, and v are all 
fixed, we can write a simple linear conversion 



Re = 0.591w - 2 - 8 



(3) 



for the range of lo we study. These low values of Re cor- 
respond to a laminar flow where both viscosity and fluid 
inertia play a role in the fluid flow. Next, to quantify the 
importance of viscous forces for influencing the horizon- 
tal motions of the beads, we consider the typical viscous 
damping time r„ ~ L 2 /v. Here we use L — l c , to quan- 
tify interactions across the length of the drum, and find 
t„ ~ 80 s, on the same order as the particle interaction 
time T co i ~ 100 s. 



III. RESULTS 



A. Trajectories 



In Ref. 0] , Mullin et al. describe three types of behav- 
iors for the three bead case. At low Reynolds number 
(Re < 1.21), they observed fixed-point behavior, where 
the beads were completely independent of each other. At 
1.21 < Rc < 2.12, the beads underwent cascading mo- 
tion in the y — z plane with the x positions fixed, with 
the axes defined as drawn in Fig. [3] It was noted that as 
they cascaded, the outer beads were stably out of phase 
with each other, and the middle bead was at an interme- 
diate position between the two. At Re = 2.12 there was 
a reversible transition to a chaotic regime where the par- 
ticles started to wander erratically in the horizontal (x) 
direction while still cascading in the y—z plane. At Re 
= 4.53, there was a transition to what Ref. @ described 
as solid body motion. In all cases, the motion in the y 
direction is always the simple cascading motion, and the 
motion in x is nontrivial; thus, like Ref. [5j, we will focus 
on the x motion for our analysis. 

In our experiment, we observed five distinct types of 
behavior within the cascade regime described by Ref. [f| , 
three of which are new behaviors. At the lowest rotation 
rates, the three beads simply cascade in the y—z plane 
with no significant motion in the horizontal (x) direction, 
corresponding to the cascading observations of Ref. [f|. 



FIG. 3: (Color online) The axes are denned such that x rep- 
resents the horizontal direction, y is the vertical, and z is the 
depth away from the front edge of the drum. Note that, with 
our current apparatus, there is no way to measure z directly. 



Figure H^a) shows a typical example of this periodic tra- 
jectory, at a rotation u) — 5.31 rad/s. Each bead appears 
to move isolated from the influence of the other beads. 
As rotation rate is increased beyond this initial simple 
periodic regime, there are several observed types of tra- 
jectories, depending on rotation rate. The two simplest 
of these are also periodic, but to differentiate their unique 
behaviors we have labeled them as doublet [Fig.^Jb)] and 
triplet [Fig.[3Jc)] states. In the original study |5J, neither 
the doublet nor triplet states were mentioned. 

In the doublet state, two of the beads will lock to- 
gether so that they are cascading in one another's wakes. 
The determination of which two beads will tend to pair 
up is a result of initial conditions, and not a systematic 
trend in the experimental apparatus. Simply stopping 
and restarting the drum at the same rotation rate can 
sometimes switch which two beads will form a pair. 

In the triplet state, all three beads come together and 
cascade in line with one another, in a similar fashion to 
the doublet state. The three beads can be stacked on top 
of each other, touching, or they can be spaced out within 
the drum, following each others' wakes without touching, 
depending on whether the beads are closer to in phase 
or out of phase as they approach one another. Both the 
doublet and triplet states are stable configurations and 
have been tested to remain locked for periods exceeding 
twenty-four hours in duration. 

At certain rotation rates, the beads will wander chaot- 
ically in the x direction. For chaotic trajectories with a 
low enough lu, there will be a bias to one side of the drum 
or the other. This biased chaotic trajectory is illustrated 
by a typical example as shown in Fig. (Ud). Two of the 
three beads will tend to repeatedly approach and interact 
with one another, while the third bead will remain seg- 
regated at the far end of the drum. This third bead still 
experiences long range forces from the other two beads, 
and can be seen to move in phase with the collisions of 
the other two. The determination of which two beads 
will tend to pair up is again seemingly a result of initial 
conditions, analogous to the doublet state. Like the dou- 
blet and triplet states, this biased chaotic behavior was 
not seen previously [5|]. 
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FIG. 4: (Color online) Typical examples of trajectories cor- 
responding to the five regimes we observe, (a) Periodic tra- 
jectory with lu — 5.31rad/s. (b) Stable doublet state with 
lu — 7.14 rad/s. (c) Stable triplet state with lu — 7.88 rad/s 
(d) Biased chaotic state with lu = 6.56 rad/s. In such states, 
the center bead strongly interacts with only one of the outer 
beads, (e) Fully chaotic state with lu — 8.51 rad/s. For all 
figures, the varying thickness of the lines is a result of par- 
ticle position uncertainty due to optical distortion from the 
curved drum walls. This distortion is exaggerated due to par- 
allax near the end caps, and minimized in the center of the 
drum. 

For higher oj chaotic states, the beads explore a more 
rich set of interactions, where they wander somewhat er- 
ratically in the horizontal direction, occasionally even 
colliding with each other. The collisions observed in- 
clude pair collisions (left- middle and right- middle) as well 
as triplet-type collisions where all three beads come to- 
gether. Figure |4je) shows a typical example of this be- 
havior, which we call the fully chaotic state, similar to 
the behaviors illustrated in Ref. 5|. In pair collisions, 
beads can be either in phase or out of phase with one 
another (in the cascading direction). In phase collisions 
are more direct, with the beads immediately colliding 
and moving away, while out of phase collisions often in- 
volve the beads cascading over one another several times 
before colliding. Triplet-type collisions generally involve 
two beads cascading over one another while a third bead 
approaches and collides with them. Note that when we 
have out of phase collisions, the two beads have identi- 
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cal x positions for a while; our data acquisition rate is 
not fast enough to carefully follow their motion in y, and 
we cannot distinguish between the beads at that point. 
Thus, it is likely that in some cases, the two beads ex- 
change places, but we cannot detect this. For example, 
in Fig. |4] the middle bead is always drawn with the same 
color, but it is important to recognize that it is quite pos- 
sible that the identity of this bead changes at collisions. 



B. Phase Diagram 

To probe the dependence of the particles' behavior on 
rotation rate, we recorded sixty-eight videos at rotation 
rates ranging from 5.1 to 12.3 radians per second. An 
analysis of these videos allows us to map out a phase di- 
agram as shown in Fig. [5] The colored blocks denote dif- 
ferent regimes, and gray blocks represent regimes where 
there is some overlap of behavior, or transitions between 
two regimes. The width of these transition blocks is due 
to uncertainty both due to the measurements themselves, 
and also to the discrete, digital motor control circuit, 
which limits the resolution of u to 1% as noted in Sec.HTl 
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FIG. 5: (Color online) Phase diagram showing the behaviors 
observed at various rotation rates. The symbols below the 
bar indicate individual observations. 

With increasing rotation rate u>, the system undergoes 
a phase transition from periodic to biased chaotic behav- 
ior at to ~ 5.7, which corresponds to a Reynolds number 
Re ~ 3.4, comparable to the transition from periodic to 
cascading motion seen in Mullin et al.'s experiment at 
Rc = 2.12 [5[. The difference in Re is perhaps due to our 
different Galilei numbers; because their viscosity was 1.2 
times larger than ours, their values for Ga are smaller by 
that same ratio. 

In our observations, the biased chaotic regime is fol- 
lowed by a long doublet regime. After this doublet 



regime, we find a small window of triplet behavior around 
lu ~ 7.8 rad/s, which begins a mixed region of behavior, 
consisting of slices of both chaotic and triplet behavior, 
extending until uj ~ 10.4. For rotation rates higher than 
those at which we find triplet behavior, we find reliable 
doublet trajectories. For high enough rotation rates, we 
should transition into the motion Mullin et al. described 
as solid-body [5j. We do not probe this regime due to 
limitations of our motor driving the rotating drum. 

This phase diagram illustrates a rich landscape of in- 
teresting regimes of particle behavior with a new level 
of detail. Specifically, Ref. [f| identified only one sim- 
ple, contiguous block of chaotic behavior, while we have 
identified multiple windows of periodic behavior embed- 
ded within large chaotic regimes, as well as previously 
unidentified periodic behaviors. Furthermore, the dis- 
tinction between two different types of chaotic behavior 
illustrates the complexity of the system. 

There are two inter-related caveats to be considered 
when discussing the phase diagram in Fig. transient 
behavior and motor drift. Transients can pose potential 
issues in situations where the transients last longer than 
the duration of an experiment. In a given experiment, 
after the drum begins rotating, the system takes some 
time to settle into its long-term behavior. For example, 
in Fig.[4jc), the particles move back and forth across the 
drum, colliding several times before finally coming to- 
gether to form the triplet state at T ~ 12 min. In this 
case, the time is small compared to typical experimental 
durations (T ~ 300 min). However, in other experiments, 
such as that shown in Fig. [51 transient behavior can per- 
sist for longer periods of time. Here, the trajectory is 
seemingly chaotic for ~ 120 min before settling into a 
stable triplet configuration. 
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FIG. 6: (Color online) In this experiment, at uo = 8.05 rad/s, 
the particle follows a seemingly chaotic trajectory for ~ 120 
min, but then settles into a stable triplet state. 

To further explore the impact of these long transient 
trajectories, we examine each trajectory and manually 
determine an approximate transient duration. Figure [7| 
shows this transient duration plotted versus the rotation 
rate of the drum. The symbols in the graph represent 
the type of trajectory found after the transient behavior 
has died out. 

Many trajectories have relatively short transient times, 
with transients rarely exceeding sixty minutes in dura- 
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FIG. 7: (Color online) The approximate time taken for the 
initial transient behavior to die out is obtained by eye and 
plotted versus rotation rate. Symbols represent the long-term 
phase of each trajectory, after its transients have died out. 



tion. However, there are also trajectories which con- 
tain much longer transient durations, with most of these 
long-duration transient trajectories clustered around the 
transitions between different phases. One possible ex- 
planation for the long-lived transients is the drum rota- 
tion rate, which as Sec. Ill Al noted is stable to within 
~ 1%. If we consider the transition around u> = 8 
rad/s, we see that, for a given trajectory, u> could vary 
from 7.92 to 8.08 rad/s. Thus a possible source for the 
long transient behaviors is drift in rotation rate of the 
drum motor. If we imagine small windows of triplet 
behavior within a chaotic regime, a drum rotation rate 
which starts within the chaotic regime could drift into the 
triplet regime, leading to a trajectory which eventually 
"finds" the triplet state. As already discussed, the triplet 
state is very robust and stable, and thus once a trajec- 
tory finds this state, it would be very difficult to break 
out of it, even if the rotation rate wanders subsequently. 
The fact that long transients tend to cluster around the 
transitions between regimes supports this hypothesis. 

This issue of motor drift also has the potential to ob- 
scure some detail in the phase diagram. The motor con- 
trol has finite resolution in available rotation rates, and 
so there may be small windows of behavior which we are 
unable to locate. Similarly, even if we did sample these 
windows, motor drift could take the rotation rate out 
of a window if it existed within a very narrow range of 
rotation rates. 



C. Qualitative Fluid Behavior 

To qualitatively describe the fluid flow within the 
drum, we added a small quantity of Kalliroscope rheo- 
logical fluid to the glycerol within the drum. Kalliro- 
scope is a water-based suspension of microscopic crys- 
talline platelets. When placed within a moving fluid, the 
platelets tend to align such that their long axis is paral- 




FIG. 8: Kalliroscope images show the fluid behavior in var- 
ious phases: (a) periodic (u) = 5.31 rad/s), (b) doublet 
(u> = 7.14 rad/s), (c) chaotic (ui = 8.51 rad/s). 



lei to the plane of shear. Thus, the platelets will reflect 
different amounts of ambient light depending on the lo- 
cal flow of the fluid. This allows us to visualize the flow 
behavior in each of the states. 

Within the periodic regime, the three beads each have 
a well defined wake, which is bounded on each side by 
swirling, vortex-like behavior rotating about an axis that 
extends in the radial direction, as seen in Fig. Ufa). At 
the midpoints between each pair of particles, there are 
well defined shear planes which span the entire height of 
the drum. 

In the doublet regime, the two beads which are paired 
up form a wake which keeps them aligned with one an- 
other. This wake is bounded on each side by vortex-like 
regions where there is swirling fluid flow, as shown in 
Fig. |SJb). The single bead, well-separated near the far 
end of the drum, also has a well defined wake, but there 
is significantly less vortex-like behavior in the fluid. 

In the triplet regime , there is one strong wake in which 
all three beads cascade (not shown). There is a large 
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amount of vortex- like swirling that bounds this wake and 
likely leads to the observed stability of the triplet state. 

The flow within the biased chaotic regime (not shown) 
appears similar to that within the periodic regime, ex- 
cept when the beads collide. As the beads approach a 
collision, their wakes overlap and partially merge. At the 
same time, as the beads are approaching one another, the 
shear plane that separates them oscillates with greater 
and greater amplitude, until it breaks up as they ap- 
proach. After a collision, when the beads are moving 
apart, the fluid to the outside undergoes a strong vortex- 
like swirling until the beads are well separated. 

Within the fully chaotic regime, the beads' wakes are 
often less well defined and more difficult to identify, 
with large regions of complicated fluid flow, as shown in 
Fig. [SJc). However, when the beads are well separated, 
their wakes are evident, with the wakes becoming mixed 
and obscured as the beads approach one another. The 
well defined shear planes seen separating the beads in 
previous cases are not evident in the fully chaotic regime. 

In all cases the visualization makes it clear that there 
is no turbulence, in agreement with the low Reynolds 
number (Re - 2 - 8, see Sec. IHCj) . 
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D. Timescales 

From Fig. HJe), we note that the particles spend the 
bulk of their time in a well-separated state where the 
three particles are spaced far apart in the drum. This 
configuration is similar to the stable configuration seen 
in the periodic state shown in Fig. Ufa) . Disturbances of 
the trajectories away from this well-separated configura- 
tion are relatively short by comparison. An interesting 
question, then, is how much time the particles spend in 
this well separated configuration. 

A visual inspection of representative chaotic trajecto- 
ries seems to indicate a typical time between collisions of 
the particles. For example, in Fig.QJd) many of the colli- 
sions between particles occur roughly 1-2 minutes apart. 
Fourier spectra of the x trajectories are noisy and do not 
depend in any obvious way on the drum rotation rate 
to. These Fourier spectra give an indication of typical 
collision timescales, with typical peak frequencies / be- 
tween 0.3 — 1.5 cycles per minute, with large changes in 
/ at nearly the same u>, corresponding to collision times 
between 0.6 and 3.3 minutes. 

To further examine the collision times, we analyzed the 
time between particle interaction events. To define these 
events, we note that deviations from the well-separated 
configuration can be identified by simply looking for local 
minima and maxima in the trajectory x-i of the middle 
bead. Figure [HI shows the distribution of the time r be- 
tween sequential peaks in xi for three representative ex- 
periments within the chaotic regime. There is no simple 
trend in the shape of these graphs with varying to. Over- 
all, the distributions of times are all broad with standard 
deviations comparable to their means, reflecting that the 



FIG. 9: (Color online) Distributions of time between "colli- 
sions" where two or three beads come close together. These 
all correspond to fully chaotic states, with the drum rotation 
rates as indicated. 



particle trajectories are unpredictable. That is, particles 
spend a significant time in well-separated positions, and 
then begin to come together for a collision event after a 
variable amount of time. Investigation of sequential pairs 
of dwell times, Tj, showed no structure, further im- 
plying unpredictability. The state at lo = 8.45 rad/s 
shows a bimodal distribution, which we observed in only 
two out of the twenty chaotic states; the other state was 
at lo = 8.91 rad/s, with several uni-modal distributions 
observed at intermediate values of lo. 



E. Reduced Dimensionality 

To this point, all analysis has focused on the horizontal 
(x) direction trajectories, and neglected cascading in the 
y — z plane. To further simplify the number of variables 
used in the data analysis, we sought a reduced dimension- 
ality set of variables which still contains the interesting 
behavior of the system. We note that, once trajectories 
have settled into their long-term behaviors, and the tran- 
sients have died out, the center of mass of the system is 
constant within the combined noise in the three particle 
trajectories. This implies that the absolute positions of 
the particles are not needed to capture the interesting 
behavior of the system, and we can use a reduced dimen- 
sionality to study the behavior. Specifically we use the 
distances of each outer bead from the center bead: 
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x 2i = \ x 2 ~ xi\/2r b (5) 
x 32 = \ x 3 - x 2 \/2r b . (6) 

Figure [TO] shows a sample trajectory comparing the orig- 
inal coordinates with the reduced coordinates. 

l5 i wi 




t (min) 



FIG. 10: (Color online) (a) Chaotic particle trajectories for 
an experiment with ui = 8.79 rad/s. Recall that we cannot 
distinguish the spheres from one another, so that by definition 
X3 > x% > x\. (b) The same trajectories in reduced coordi- 
nates. The minima in these simplified trajectories indicate 
points where particle pairs approach one another. 

As shown in Fig. 1101 there are certain configurations 
in which the system spends more time, and other con- 
figurations that are only visited briefly. Configurations 
with the three beads well separated appear to be most 
common, while situations where the beads are close to- 
gether are shorter-lived. To quantify this, we plot a two 
dimensional histogram of the configurations, which gives 
us a way to visualize the relative amount of time each 
particle spends in various regions of phase space. His- 
tograms are plotted with a logarithmic intensity, so rare 
events can still be identified. In each experiment, we vi- 
sually inspect the data set and remove any obvious initial 
transient behavior manually before analysis. For exam- 
ple, the analysis of a triplet data set only includes the 
time after the three beads have lined up. 

Figure [TT1 shows an example of one of these histograms 
for the same experiment shown in Fig.[TOJ Notice that the 
darkest red region is in the area around x' 21 ~ x' Z2 ~ 5 — 6, 
corresponding to a configuration where the three beads 
are spread far apart, and spaced roughly equidistantly. 
There are also small clusters at x' 21 ~ 7.5 and x' 32 ~ 0, 
and its mirror x' 21 ~ and x' 32 ~ 7.5, which correspond 
to configurations where two beads are close together, and 
the third bead is far away. Finally, there is another faint 
cluster at x' 21 ~ and x' 32 ~ 0, corresponding to a state 
where all three beads are grouped together. The faintncss 
of this cluster implies that very little time is spent in this 
configuration. 




2 4 6 8 
x n 

FIG. 11: (Color online) A typical two dimensional histogram 
for a chaotic trajectory at uo — 8.79 rad/s. The axes represent 
the distances between the pairs of particles, normalized by the 
particle diameter, and color represents the number of points 
that were counted in each bin. Darker colors indicate more 
prevalent configurations. This experiment has a calculated 
entropy S = 6.33. 



The histograms show how the different phases of the 
system explore phase space. Figure [T2T a) shows a 
schematic of typical histograms for periodic, doublet, 
and triplet trajectories. Periodic trajectories have well- 
separated beads and thus the weight of the histogram 
stays concentrated at point 1. Doublet states have a 
histogram with points clustered tightly near one of the 
locations marked 2; for example, if beads 2 and 3 are 
together in the doublet state, then x' 23 = 0. In the 
triplet state, the three beads coincide in their x coor- 
dinates and thus the histogram is at the origin, where 
point 3 is shown. Figure [T2Tb) shows a biased chaotic 
trajectory, where the beads spend time in the same re- 
gion of phase space as the periodic state, but also wander 
chaotically in the x direction, smearing the histogram in 
the direction of the bias (in this case, beads 2 and 3 come 
close together). In Fig. [T27 c) we see a chaotic trajectory. 
The beads explore the phase space around the locations 
of periodic, doublet, and triplet states [points 1, 2, and 
3 in Fig. 112( a)]. Time is also spent in an intermediate 
configuration (x' 12 ~ x' 23 « 3.5) where the three beads 
are closer together than in the periodic trajectory, but 
not clustered like in the triplet state. This is in contrast 
to the chaotic state shown in Fig.[T2Td). where the beads 
spend the bulk of their time well-separated, as in the pe- 
riodic state, with occasional doublet-like collisions. The 
fact that these states occur at rotation rates which are 
very similar (lo = 9.25,9.36 rad/s) illustrates how sensi- 
tive the experiment is to the rotation rate. 
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F. Entropy 

To quantitatively study the extent to which a given tra- 
jectory explores phase space, we define a configurational 
entropy based on these histograms. If we first normalize 
a given histogram so that the sum of all bin values is 
equal to unity, the histogram will represent a probability 
distribution P with matrix elements Pij . We then define 
the entropy 

S = -k Py lnPij , with k = 1. (7) 

ij 

This entropy value is calculated for each individual tra- 
jectory and plotted versus rotation rate in Fig. []2J The 
triplet states have the lowest entropy, followed by the 
periodic, doublet, biased chaotic, and finally the fully 
chaotic states with the highest entropy. To quantify er- 
ror in entropy measurements we split each trajectory into 
two halves, and calculate the entropies Si and S2 for each 
half independently. Then the error can by defined as: 

<rs = \S 1 -S 2 \, (8) 

yielding the error bars shown in Fig. 1131 




02468 02468 

FIG. 12: (Color online) Histograms for each phase of behav- 
ior clearly illustrate the amount of phase space they explore, 
(a) Schematic of typical areas occupied by the the non-chaotic 
trajectories. "Periodic" trajectories such as shown in Fig.|4ja) 
appear as a tight cluster of points near position 1. "Doublet" 
trajectories reside at one of the two locations marked as po- 
sition 2. "Triplet" trajectories have all the particles at the 
same x position, and thus the histogram for these states is 
a tight collection of points near the origin, marked as posi- 
tion 3. (b) A histogram for a biased chaotic trajectory at 
uj — 6.05 rad/s. For this state, the entropy S = 3.52. (c,d) 
Histograms for chaotic trajectories at uj — 9.25,9.36 rad/s, 
respectively. The entropies of these states are S = 6.82, 6.49. 
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FIG. 13: (Color online) The entropy for each rotation rate, 
calculated from the 2-D histograms. 



As could be expected, the entropy is lowest in the sim- 
plest states: the triplet, doublet, and periodic states. 
Nonzero values for the entropy correspond to slight wob- 
bling of the particles around their mean positions in each 
state. At values of uj corresponding to transitions be- 
tween states, Fig. [T3] shows sharp changes in the entropy. 
For example, the biased chaotic state has entropy values 
markedly higher than the adjacent doublet states. The 
fully chaotic states have the highest entropy. Within each 
regime, there can be moderate fluctuations of the entropy 
values, with little systematic dependence on uj. 

Much as the histograms provide a visual indication of 
the degree to which a trajectory explores phase space, 
these entropy values provide a qualitative measure of that 
exploration. Low entropy periodic, triplet, and doublet 
regimes do not explore phase space much, as expected. 
Chaotic behaviors, on the other hand, have large en- 
tropies, corresponding to extensive exploration of phase 
space. These entropies vary little by comparison to the 
large jumps between regimes, indicating that the amount 
of phase space explored by chaotic behaviors does not 
vary much with rotation rate. 

The magnitude of these entropy values varies with the 
size of bins used for the histograms. We divide each axis 
into 60 divisions, giving a total of 60 2 bins, with each 
bin having a width w ~ rb/3. The qualitative results are 
unchanged with 40 2 or 80 2 bins. 



IV. SUMMARY 

We have studied a geometrically simple system con- 
taining three particles moving within a fluid-filled rotat- 
ing drum which yields a rich and varied set of behaviors. 
The phase diagram for this system showed five types of 
behavior. The first is a periodic regime where the beads 
simply cascade in the y — z plane. The second is a pre- 
viously unreported biased chaotic regime where two of 
the beads wander chaotically in the horizontal x direc- 
tion and collide with one another. The third is a doublet 
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regime, where two beads pair up and cascade on top of 
one another while leaving the third bead behind. There 
is also is a mixed chaotic regime spanning a wide range 
of rotation rates, where the beads wander chaotically in 
the horizontal direction, with all three beads interacting 
and mixing. Within this mixed regime, there are small 
windows of rotation rates which result in triplet behav- 
ior, where all three beads will line up and cascade on top 
of one another. Finally, we find a regime where triplet 
behavior is the only type of trajectory seen. 

The question of transient behavior deserves significant 
attention and exploration. We find chaotic states that 
are persistent over many hours. However, based on sim- 
ilarities between the chaotic states and the transient be- 
havior of the triplet states, it is possible that the chaotic 
states could eventually fall into a stable triplet state. 
Our "fastest" observed chaotic states have mean collision 
times of approximately 0.6 min, and our longest observa- 
tions are up to 360 min, so at most we have 600 collisions 
observed in chaotic states without a transition to a triplet 
state, thus at least showing that if these are transients 
that they are very long-lived. 

Another interesting question is that of dependence on 
initial conditions. Due to the way the particles were po- 
sitioned within the drum, it was difficult to control their 



exact starting positions. However, the long-term statis- 
tics of the system's behavior were reproducible, within 
the limits of motor drift and transients. 

We have also tried preliminary experiments with a 
longer drum (thus larger aspect ratio) , and find that par- 
ticles prefer doublet or triplet states; we do not see long- 
lived chaotic states in longer drums. This suggests that 
a key to the chaotic states is indeed the finite drum size 
that forces the beads to interact with each other. Overall, 
we have demonstrated a simple system with complex and 
chaotic behavior, and furthermore shown that within the 
previously observed chaotic regime (Ref. @) there is a 
richness of behavior with the "degree" of chaos changing 
dramatically with only slight changes of rotation rate. 



Acknowledgments 

We thank M. Schatz and D. Borrero for guidance in 
the use of Kalliroscope and many helpful discussions, and 
T. Mullin for his inspiration of this project and helpful 
discussions. This work was supported by the Emory Uni- 
versity Graduate School of Arts and Sciences and NSF 
DMR-0804174. 



H. Shin & M. R. Maxey. "Chaotic motion of nonspherical [11 
particles settling in a cellular flow field." Phys. Rev. E, 
56, 5341-5443 (1997). 

Anil, Satheesh, & T. R. Ramamohan. "Chaotic dynamics [12 
of periodically forced spheroids in simple shear flow with 
potential application to particle sedimentation." Rheo- 
logical Acta, 34, 504-511 (1995). 

H. Aref & S. Balachandar. "Chaotic advection in a stokes [13 
flow." Phys. Fluids, 29, 3515-3521 (1986). 

K. Asokan, Anil, J. Dasan, K. Radhakrishnan, Satheesh, 
& T. R. Ramamohan. "Review of chaos in the dynam- [14 
ics and rheology of suspensions of orientable particles in 
simple shear flow subject to an external periodic force." 
J. Non-Newtonian Fluid Mech., 129, 128-142 (2005). [15 
T. Mullin, Y. Li, D. C. Pino, & J. Ashmore. "An ex- 
perimental study of fixed points and chaos in the motion [16 
of spheres in a stokes flow." IMA Journal of Applied 
Mathematics, 70, 666-676 (2005). [17 
D. D'Humieres, M. R. Beasley, B. A. Huberman, & 
A. Libchaber. "Chaotic states and routes to chaos in 
the forced pendulum." Phys. Rev. A, 26, 3483 (1982). [18 
S. J. Van Hook & M. F. Schatz. "Simple demonstrations 
of pattern formation." The Physics Teacher, 35, 391-395 
(1997). [19 
P. N. Segre, E. Herbolzheimer, & P. M. Chaikin. "Long- 
range correlations in sedimentation." Phys. Rev. Lett., 
79, 2574 (1997). 
K. O. L. F. Jayaweera, B. J. Mason, & G. W. Slack. [20 
"Behavior of small clusters of spheres falling in a viscous 
fluid." J. Fluid Mech, 20, 121 (1964). [21 

I. M. Janosi, T. Tel, D. E. Wolf, & J. A. Gallas. "Chaotic 
particle dynamics in viscous flows: The three-particle 
stokeslet problem." Phys. Rev. E, 56, 2858 (1997). 



M. Tirumkudulu, A. Tripathi, & A. Acrivos. "Particle 
segregation in monodisperse sheared suspensions." Phys. 
Fluids, 11, 507-509 (1999). 

M. Tirumkudulu, A. Tripathi, & A. Acrivos. "Erratum: 
"particle segregation in monodisperse sheared suspen- 
sions" [phys. fluids [bold 11], 507 (1999)]." Phys. Fluids, 
11, 1962 (1999). 

K. Choo, T. C. A. Molteno, & S. W. Morris. "Travel- 
ing granular segregation patterns in a long drum mixer." 
Phys. Rev. Lett, 79, 2975-2978 (1997). 
K. Choo, M. W. Baker, T. C. A. Molteno, & S. W. Mor- 
ris. "Dynamics of granular segregation patterns in a long 
drum mixer." Phys. Rev. E, 58, 6115-6123 (1998). 
J. Williams. "The segregation of particulate materials, a 
review." Powder Technology, 15, 245-251 (1976). 
J. Bridgwater. "Fundamental powder mixing mecha- 
nisms." Powder Technology, 15, 215-236 (1976). 
K. M. Hill & J. Kakalios. "Reversible axial segregation 
of binary mixtures of granular materials." Phys. Rev. E, 
49, R3610-R3613 (1994). 

K. M. Hill, A. Caprihan, & J. Kakalios. "Bulk segrega- 
tion in rotated granular material measured by magnetic 
resonance imaging." Phys. Rev. Lett., 78, 50-53 (1997). 
P. N. Shankar & M. Kumar. "Experimental determi- 
nation of the kinematic viscosity of glycerol-water mix- 
tures." Royal Society of London Proceedings Series A, 
444, 573-581 (1994). 

E. R. Benton & A. Clark. "Spin-up." Annual Review of 
Fluid Mechanics, 8061, 257-279 (1974). 
H. D. Baehr & K. Stephan. Heat and Mass Transfer 
(Springer) (2006). 



